Spontaneous-emission rates in finite photonic crystals of plane scatterers 
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The concept of a plane scatterer that was developed earlier for scalar waves is generalized so 
that polarization of light is included. Starting from a Lippmann-Schwinger formalism for vector 
waves, we show that the Green function has to be regularized before T-matrices can be defined in a 
consistent way. After the regularization, optical modes and Green functions are determined exactly 
for finite structures built up of an arbitrary number of parallel planes, at arbitrary positions, and 
where each plane can have different optical properties. The model is applied to the special case 
of finite crystals consisting of regularly spaced identical planes, where analytical methods can be 
taken further and only light numerical tasks remain. The formalism is used to calculate position- 
and orientation-dependent spontaneous-emission rates inside and near the finite photonic crystals. 
The results show that emission rates and reflection properties can differ strongly for scalar and for 
vector waves. The finite size of the crystal infiuences the emission rates. For parallel dipoles close to 
a plane, emission into guided modes gives rise to a peak in the frequency-dependent emission rate. 

PACS numbers: 42.70.Qs, 02.30.Rz 78.66.-w, 42.50.-p 



I. INTRODUCTION 

Photonic crystals are a well-studied subject nowadays, 
both theoretically and experimentally 0]. Of funda- 
mental importance is the prediction that in three- 
dimensional photonic crystals that meet a tough combi- 
nation of requirements, light propagation will be com- 
pletely inhibited in all directions and a photonic band 
gap will show up for certain frequencies of light. It is 
important for technology that photonic crystals can be 
created that guide light with low losses, and bend light 
on a scale of an optical wavelength. The latter properties 
do not require a band gap in all three dimensions. 

A photonic-band-gap crystal would reflect light for all 
angles of incidence, when the frequency of the light lies 
within the gap. However, lower-dimensional photonic 
crystals such as Bragg mirrors can also be omnidirec- 
tional mirrors, without having a band gap Thus, 
external light sources can only give an indication that 
there is a band gap or a proof that there is no gap. 

Internal light sources such as excited atoms do a better 
job in probing a band gap, because only a gap would com- 
pletely inhibit spontaneous emission by internal sources 
For the same reason, a photonic-band-gap crystal 
would be a whole new playground in quantum optics, 
both when one is interested in spontaneous emission in 
itself, and in processes which normally are obscured or 
made less efficient because of spontaneous emission. Not 
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only emission rates would be strongly modified inside 
a band-gap crystal, but also resonant dipole-dipole in- 
teractions, for example, as they are mediated by the 
electromagnetic field 0. The focus of this paper is on 
spontaneous-emission rates of visible light. 

For atomic transition frequencies in the band gap of 
an infinite three-dimensional photonic crystal, emission 
rates vanish everywhere in the inhomogencous structure. 
In practice, such a uniform suppression of emission rates 
has not yet been observed for visible light: evidence of 
crystals exhibiting a full photonic band gap in the visible 
has not been reported to date. Even when in the fu- 
ture such crystals will exist, position-dependent emission 
rates will occur at the edges of the crystals. In general, 
spontaneous-emission rates of inhomogencous dielectrics 
with a high refractive- index contrast, including photonic 
crystals, are strongly position and orientation dependent. 
Calculated spontaneous-emission rates in this paper will 
prove this point. Also, finite-size effects will show up in 
our calculations. The model studied here is a finite pho- 
tonic crystal consisting of a finite number of parallel and 
infinitely thin planes. More about our model will be said 
later in this Introduction. 

In many experiments, dipole orientations are hard 
to control. When averaged over dipole orientations, 
spontaneous-emission rates are proportional to a quan- 
tity called the 'local optical density of states' (LDOS) 
0, ■ The concept of a local density of states was bor- 
rowed from solid-state physics. The local optical density 
of states was first named the 'local radiative density of 
states' 0, which is the same quantity. 

Interestingly, the calculation of position-dependent 
spontaneous-emission rates also has a bearing on the in- 
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terpretation of measurements performed with a near-field 
scanning optical microscope (or NSOM). In these mea- 
surements, a sample is illuminated through the tip of the 
microscope, and scattered light is recorded. In a simple 
model, the disturbance of the optical field by bringing 
the tip of the microscope to the sample is assumed to 
be weak, and the tip is modelled as a dipole with a cer- 
tain strength and orientation. Then, if the light scattered 
in all directions would be recorded, the measured signal 
would be proportional to the local spontaneous-emission 
rate at the position of the tip in the absence of the tip 
i 

After the above general considerations, we now turn 
to the topic how spontaneous-emission rates inside pho- 
tonic crystals are actually being calculated. The exis- 
tence of a band gap in an infinite photonic crystal can 
be inferred from a band structure calculation, which for 
a three-dimensional photonic crystals is an art in itself 
(see the recent review ilDj). Quite another and more 
difficult matter is it to calculate emission rates inside in- 
finite crystals prlll^. Emission rates inside or near finite 
photonic crystals are even harder to calculate. Other in- 
teresting quantities would be near-field or far-field spec- 
tra of internal sources, or dipole-dipole interactions and 
superradiance effects of atoms embedded in a finite three- 
dimensional photonic crystal, to name a few complex pro- 
cesses in a complex environment. In such cases, results 
of calculations are hard to check and - even if correct - 
they might not give much insight. 

It is therefore very useful to study complex processes in 
simplified models for photonic crystals. Widely used is 
the so-called quasi one-dimensional model (or isotropic 
model) for photonic crystals where it is assumed 
that the red edge of the stop bands of the crystal occur at 
the same band-edge frequency for all three-dimensional 
propagation directions, and similarly for the blue band 
edge. Such a model will describe qualitatively correct the 
processes well inside the band gap, while overestimating 
effects of the photonic crystal at the edges of the gap. 
The isotropic model also neglects all position and orien- 
tation dependence of emission rates outside the band gap. 
Inspired by the model calculations, more realistic numer- 
ical calculations have recently appeared that indeed show 
the weaknesses of the isotropic model ^2 ■ 

In this paper another simple model is proposed, one 
which takes into account the strong spatial and orien- 
tational dependence of optical properties and the finite 
size of the crystals. On the other hand, it gives up the 
existence of a full band gap, as only variations of the 
refractive index in one dimension are considered. Dielec- 
tric slabs are modelled as infinitely thin planes, which 
will be called plane scatterers. A multiple-scattering for- 
malism is set up in which optical modes and the Green 
function (a tensor, really) can be calculated exactly for 
crystals consisting of an arbitrary number of plane scat- 
terers. The present model is a generalization of previous 
work that treated scalar waves only Q . The inclusion of 
polarization of light will turn out not to be straightfor- 



ward. 

Infinitely thin planes were used as model systems in 
photonics before, for example in |l4j where light propaga- 
tion was considered in one dimension only and for an infi- 
nite crystal. The model was generalized in (ISi] . where in- 
finite photonic crystals were built of infinitely thin planes 
and their band structure was determined for waves prop- 
agating in three dimensions. In 16], both infinite and 
finite crystals of planes were considered and their trans- 
mission and reflection properties were studied with the 
use of transfer matrix methods. The infinite crystal was 
again considered in 0,0,0,113 and named the 'Dirac- 
comb superlattice'. Frequency-dependent emission rates 
were determined for several positions in the unit cell and 
both TE- [Hill and TM-waves [Ull^ were considered. 
The periodicity of infinite crystals makes that Bloch's 
theorem can be used in the analysis. Finite photonic 
crystals do not have this advantage and the analysis of 
the optical properties is usually more difficult. Position- 
dependent spontaneous-emission rates remain to be ex- 
plored in model photonic crystals consisting of a finite 
number of plane scatterers, where light propagation in 
all three dimensions is taken into account, for both po- 
larization directions. 

It has been known for a long time that spontaneous- 
emission rates of atoms change when positioned at dis- 
tances on the order of the wavelength of light away from 
a mirror |^ |23 |. A recent surprise was the mea- 
surement and analysis that even a distant mirror (25cm 
away) can change emission rates when lenses are used 
and when the atomic positions are controlled with a sub- 
wavelength precision (a few nanometers) p^ . l25l | . 

More complicated than emission near a mirror is it 
to calculate emission rates of atoms in or near one- 
dimensional photonic crystals. A multi-purpose formal- 
ism for calculating optical modes in layered dielectrics 
|26l | was used in jj] to calculate emission rates inside fi- 
nite periodic layered structures, especially inside struc- 
tures that reflect light incoming from all directions, the 
so-called 'omnidirectional mirrors' 3]. 

Interestingly, light transmission through finite one- 
dimensional photonic crystals can be found exactly in 
terms of the transmission through a unit cell, the number 
N of unit cells, and in terms of the Bloch wave vector of 
the corresponding infinite crystal structure. In p^ . this 
is shown for a simple unit cell containing two layers, but 
it was also proven for general unit cells |23| • This remark- 
able result was reviewed in [l^, where its importance is 
stressed not only in optics but also in acoustics, quan- 
tum mechanics, and other branches of physics. In the 
T-matrix formalism of this paper (which differs from the 
usual transfer matrix method for layered dielectrics), we 
find similar analytical results, also involving the Bloch 
wave vector. Such attractive analytical results are not 
available for more complex dielectrics such as finite two- 
[30I I3H or three-dimensional photonic crystals [s^, so 
that in those cases the use of efficient numerical tech- 
niques is essential. 
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The advantage of our plane-scatterer model is that 
modes and Green functions (and therefore emission rates) 
can be calculated exactly in a Lippmann-Schwinger for- 
malism, for every finite crystal size, and that light propa- 
gation in all directions is taken into account. Lippmann- 
Schwinger formalisms arc more commonly used |3,1 | . but 
when the finite volumes of dielectric scatterers are fully 
taken into account, numerical discretization of the dielec- 
tric is required and the model stops being simple |3lll3^. 
To be sure, the simplicity of our model entails that in 
some aspects it becomes less realistic, as will be stressed 
where appropriate. 

In Sec. m multiple scattering of light is introduced 
and central equations are derived in representation- 
independent notation. In Sec. IIIII the free-space Green 
tensor is regularized and a T-matrix of a plane scatterer 
for light waves is derived. Sec. IIVI discusses all optical 
modes (propagating and guided modes, including polar- 
ization) that exist in crystals of plane scatterers. Po- 
sition and orientation dependent spontaneous-emission 
rates are calculated in Sec. [V| Conclusions can be found 
in Sec. |Vl| 

II. MULTIPLE-SCATTERING THEORY FOR 
VECTOR WAVES 

Some important equations of multiple-scattering the- 
ory |33| will be presented, mostly in representation- 
independent notation, for light in arbitrary inhomoge- 
neous dielectrics. In later sections a particular class of 
dielectrics will be studied and a suitable representation 
is chosen, but then the involved notation might obscure 
the general structure of the equations. 

The wave equation for the electric field Eo(r, uj) in free 
space is 

(u;/c)2Eo(r, tj) - V x V x Eo(r, u) = ^ (la) 
{(cj/c)2| • -V X V X }Eo(r, uj) = 0. (lb) 

The symbol I denotes the unit tensor in three- 
dimensional space. The solutions of Eq. lla() arc plane 
waves with wave vector k and polarization direction nor- 
mal to k. With the free-space wave equation (|la|l . a 
Green tensor (or dyadic Green function) is associated 
that satisfies 

{{ij/cf\--V xV X }Goir,r\uj) ^S^{r-r')\. (2) 

Let L(r, a;) be the quantity between curly brackets in 
Eqs. Ijlb|l and Both equations can be considered 

as the real-space representations of an abstract tensor 
operator L(tj) operating on the vector field Eo{uj) and 
on the Green function Gq{uj), respectively: 

L(w) • Eo(tj) = 0, L(w) • Go(w) = I® I. (3) 

The identity operator in real space is denoted by 1 and it 
has the property (r|l|r') = S^{r — r'); confusion with the 



unit tensor I should not arise; the (g) denotes the tensor 
product. 

In the presence of an inhomogeneous dispersive linear 
dielectric, the wave equation for the electric field is mod- 
ified into 

L(w) • E(w) ^ V(cj) • E(w), (4) 

where the frequency-dependent optical potential V is de- 
fined in terms of the dielectric function e(r,a;) as 

(r|V(Lj)|r') =-[e(r,w)-l](Lj/c)2|5(r-r'). (5) 

The delta-function on the right-hand side defines the po- 
tential as a local quantity (which the T-matrix, to be 
defined shortly, is not). In other words, this delta func- 
tion appears for any potential. 

The electric field Eo{lu) is modified into E(ijj), and the 
two fields are related through the Lippmann-Schwinger 
(LS) equation 

E(w) = Eo(w) -I- Go{uj) ■ V(w) • E(w) (6a) 
= Eo(w)-hGo(w) • V(a;) •Eo(t^) 

+Goiuj) ■ y{uj) ■ Go{uj) ■ y{uj) ■ Eo(t^) -I- . (6b) 
= Eo(cj) + Go(t^) • T(co) • Eo(tj). (6c) 

One can check that indeed the field E(aj) that satisfies 
Eq. (|6a|l is also a solution of Eq. Q. The solution of 
Eq. ^ can be found iteratively in higher and higher 
orders of the optical potential V, as given by the multiple- 
scattering series in Eq. I|6bp : the (dyadic) T-matrix in 
Eq. I|6c|) by definition is the formal sum of the infinite 
summation in Eq. Ij6b|l . The T-matrix is a 3 x 3 tensor. 
By combining Eqs. Ij6bp and H6c|l . the formal solution for 
the T-matrix is 

T{uj) = • [1 (g) I - Go(w) • y{uj)r^ . (7) 

The scattering problem is solved exactly once the T- 
matrix is known. 

There may exist optical modes that are bound to the 
scatterer. Such bound modes correspond to solutions of 
the LS equation (|6a|l in the absence of an incident field; 
with Eq. Q we can rewrite this homogeneous equation 
as 

T-\uj) -yiuj) -Eiuj) ^0. (8) 

It follows that bound solutions of the electric field will 
correspond to the poles of the T-matrix. Actually, 
Eq. also shows that a nonzero solution for E(w) can 
only occur when T(w) has a pole. The T-matrix not only 
solves the scattering problem for incident fields but also 
contains all information about bound modes. 

In the presence of the dielectric the Green function also 
changes, from Gq to G. The latter satisfies 

[L{lj)-\/{lu)]-G{uj)^1<»\, (9) 
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The solution for the Green function analogous to Eq. © 
for the electric field is the three-dimensional Dyson- 
Schwinger equation 



G{uj) = Go{lj) 



Go(t^) • V(w) • G(w) (10a) 
GoH-TH-GoM- (10b) 



It can be verified that a solution of (|10a|l also is a solution 
of Eq. 0. The problem how to find such a solution is 
solved once the T-matrix (0) is determined, because an 
iteration of Eq. (|10a|) analogous to the series expansion 
(|6b|l for the electric field shows that the Green function 
can also be expressed in terms of the T-matrix, as given 
by Eq. (Tmhll . 

Equation (|10|l also holds when the total potential V(ijj) 
is a sum of single-scatterer potentials Vq (lo) . By iterating 
one finds that the total T-matrix for an arbitrary number 
N of these scatterers is 



N 



=^V„-|-^V0.Go-V„+^ V^.Go-V^-Go-V„+.. 

(11) 

Often it is more convenient make an equivalent expansion 
in terms of the single-scatterer T-matrices |3^ : 

N 

Q=l a,l3{^a) 

+ ^ • Go ■ T/3 • Go ■ Ta ^ 

a,/3(#a),7(#/3) 



The frequency dependence was suppressed in Eqs. Ijllfl 
and (1121). The form Eq. (O of the total T-matrix will 
be used later in this paper, for model systems where the 
infinite summation can be performed explicitly. 



III. PLANE SCATTERERS FOR VECTOR 
WAVES 

The general results of multiple-scattering theory that 
were presented in Sec.^will now be applied to dielectrics 
that can be described as a collection of parallel planes. 
A suitable representation is chosen, and specific forms of 
the potential V, the free-space Green function Go, and 
the incoming electric field Eq are determined. With this, 
T-matrices for a single plane and for an arbitrary number 
of planes are derived. 



A. Dyadic Green function in plane representation 

A solution for the free-space dyadic Green function 
can be found in three-dimensional Fourier space. By 
translational invariance, (k|Go(w)|k') must be equal to 
(27r)353(k- k')Go(k,w). The Green function Go(k,w) 
satisfies 



{[{uj/cf-k^] I + fc2kk}-Go(k,w) = I. 



Here, k denotes a unit vector in the direction of the 
wave vector k. Equation H13|l is a 3 x 3 matrix equation 
whose representation diagonalizes in the polarization ba- 
sis {k, CTi, CT2} with the longitudinal direction k and two 
orthogonal transverse directions (Ti,2- The solution of 
(|T^ is 



G^o'(k,^) 



1 



G!^''(k,c^) 



(c/co) 



(14) 

where j denotes ci or a2- AH six non-diagonal elements 
of the Green tensor are zero in this representation. This 
is the retarded Green function once we assume that the 
frequency lo has an infinitesimally small positive imagi- 
nary part. 

The above Fourier representation is not what we need. 
It is convenient to work in the "plane representation" : 
in two-dimensional Fourier space in the directions par- 
allel to the planes and in real space in the z-direction 
perpendicular to the planes. For the polarization repre- 
sentation choose the orthonormal basis {sk, Vk, z}. Here, 
z is the unit vector in the z-direction; Vk is the unit vec- 
tor in the direction of the projection of the wave vec- 
tor k on the plane, so that the wave vector k has a 
component fc|| in the Vk-direction and its full represen- 
tation is (0,fc||,fcz); the Sk-polarization direction is or- 
thogonal to the optical plane that is spanned by the 
other two basis vectors. Then the operator L(a;) has 
,(12) the form (k,,, z|L(w)|k'||, z') = i2n)^5^{\n\ - k'\\)6{z - 
2;')L(k|| , z,uj), where the operator L(k|| , z, w) has the ma- 



trix representation 








(15) 
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The Green function in the same representation becomes 
(k||,z|Go(u;)|k'||,z') = [2nf5^^-\^\\)Qo{S,z,z',uo), 
and the transformed Eq. (|13|l is a system of differential 
equations: 



L(k||,z,a;) 



GZ 
n 



CZQ 















1 =<5(z-z') U 


1 
















(16) 

The Gq^ are the components of Gq and their arguments 
(k||,z,z',ix') were dropped for brevity. By choosing the 
plane representation, the matrix elements of Go only de- 
pend on the magnitude and not on the orientation of k|| . 
All components involving an s-label are zero, except the 
ss-component. Gff' satisfies the same differential equa- 
tion as the Green function 170 of the Helmholtz equation 
for scalar waves, so that for lo > we have 



Go''m\,z,z',uj) =5o(k||, z,z',w) 



rjikz \ z — z' I 



(17) 



(13) 



The variable k^ is not independent from k|| , but rather an 
abbreviation for [(w/c)^ — fcy]^/^. The remaining coupled 
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differential equations of Eq. (|16|) can also be solved (again 
for CO > 0), now that Gff is known: 



GS"(k||,2,z',C.) = 



90 



(18a) 



Gg'^(k||,z,z',c.) 

Gg^(k||,2,2',C.) 







■go sign(z - z') (18b) 
(18c) 

A:f5o + S{z - z')] (18d) 



Green functions in the right-hand sides are understood 
to have the arguments (k||,z,z',ti;). The above method 
of solving differential equations docs not give a value for 
the sign-function when z is equal to z' . The Green func- 
tion components (|18|l can alternatively be found from an 
inverse Fourier transformation 



ik'^ iz — z') 



1 f°° 

Go(k||,z,z',w) = — / dfc^ Go(k||,fc^,,w) e' 

(19) 

This integration can only be performed in a representa- 
tion that does not co-rotate with k'^ . The basis of Eq. H14|) 
is not adequate, but again the basis {sk, Vk, z} suits well. 
With this Fourier method one finds the value for the 
sign-function in Eq. Ijl8b|l when z equals z': for z — z' 
the relevant integrands in Eq. (|19l) are antisymmetric in 
the variable k'. 



B. Regularization of the Green function 

The T-matrix of a plane scatterer for vector waves can 
be found by solving the appropriate Lippmann-Schwinger 
equation (Q. A plane wave incident from z = — oo with 
wave vector k and arbitrary amplitude Eq and trans- 
verse polarization vector cTk = (cTg, (t„, cr^) is scattered 
by a plane at z = Zq. Because of the symmetry in the 
in-plane directions, it is convenient to choose the plane 
representation for the LS equation. In terms of the Dirac- 
notation, the electric field is a "ket" ; the plane represen- 
tation is found by taking the inner product of Eq. ® for 
the electric field with the "bra" (k|[, z|, and by inserting 
the unit operator 



(2^) 



L_ j dViidz' |k'||,z')(k'||,z'| 



(20) 



A plane is assumed to be infinitely thin and it can be de- 
scribed by the optical potential \/(z, u) = V{uj)5{z— Za)\- 
(A specific model potential will be chosen in Sec. lIIIFl ) 
The integral can be evaluated immediately and we get 



at the positions of the dots in the representation- 
independent equation (0. The incident field takes the 
form EkCT,o(k|| , z, w) = i?of k exp(ifc2z). The solution of 
the LS equation corresponding to this incident field is 
EkCT(w). The LS equation in the mixed representation 
becomes 

Ek^(k||,z,u;) = £;ocrke^'=^^ (21) 

/+00 
dz' Go(k||, z,z', w) • V(z', w) • Eko-(k|| , z', w). 
-OO 



EkCT(k||,z,u;) 



Eorrv^e'^'^ (22) 
V^(w)Go(k||, z,Za,w) • Eko-(k||, Zq, w). 



The usual way to solve this equation would be to put the 
position z equal to Zq, and to solve for Ek<T(k|| , Zq, cj). 
The result would then be inserted back into the above 
equation to obtain an expression for EkCT(k|| , z, w). 

However, the Green tensor Go is not defined when the 
positions z and Zq are identical, because of the delta func- 
tion in the component Gq^ [Eq. (|18d|l ]. One could just 
neglect the delta function, as might be correct in other 
situations [H^Ij but it will be argued in Sec. IIII Cl that 
this procedure would be wrong in our case. Therefore, 
the Green tensor (|18|l is not suited for setting up a the- 
ory for the scattering of vector waves by infinitely thin 
planes. 

It is known that "regularization" of Green functions is 
sometimes needed when modelling finite-sized scatterers 
as mathematical objects with zero volume, in order to 
have a model that is relevant for optics. (Not always: 
regularization was not needed for scalar waves scattering 
off planes ^\.) In a regularization procedure usually a 
cutoff parameter is introduced that modifies the behav- 
ior of Green functions at distances much smaller than op- 
tical wavelengths, and mathematical problems arc thus 
overcome. In some cases, the regularization parameter 
can be sent to infinity in the final stage, while in other 
cases the cutoff parameter must be kept finite. For exam- 
ple, for point scatterers the problem of diverging Green 
functions occurs both for scalar and for vector waves. 
Point scatterers have been studied extensively and sev- 
eral regularization schemes have been proposed (see 
and references therein). 

The same regularization procedure will now be cho- 
sen for plane scatterers as was done before for point 
scatterers j^: a high-momentum cutoff is introduced 
in three-dimensional Fourier space: instead of the free- 
space Green function Go(k, cj) of Eq. H14|l . a regularized 
free-space Green function Go(k, uj) will be used. The lat- 
ter is defined in terms of the former as 



Go(k, w) 



A2 



A2 



Go(k,cj) 



(23) 



The cutoff momentum A is assumed to be much larger 
than the magnitude oj/c of the optical momentum, so 
that at optical wavelengths Go — Go. The effect of this 
cutoff in the real-space representation is also known [s^ . 
Here its effect on the Green function in the plane repre- 
sentation is important. After an inverse Fourier transfor- 
mation, again only in the z-direction, one obtains (again 
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for ll> > 0) 

Gg"(k||,z,zi,tj 

GS^(k||,z,zi,c. 
GS^(k||,z,zi,cj 

Go"(k||,2:,Zi,C^ 

Gg^(k|j,z,zi,a; 



_ A2 / g-A|||z-^;i 

A2 + {uj/cj^ ^ 2A^ 

1,2 2 

A2sign(z — zi) k\\ 
A2 + (u;/c)2 (w/c)2 

fc,go + (*/2)e-^iil^-^il 







A2 



X 



A2 + (t^/c)2 
fc2c2 A2 + (c./c)2 ^ , , 



24a) 
(24b) 

(24c) 
(24d) 

(24e) 



2A||(c^/c)2 



In the right-hand sides, the arguments (k|| , z, zi,uj) of the 
Green functions were dropped; A|| is short-hand notation 
for (A2 + fc2)i/2. again, the sign-function is zero when its 
argument is. 

All components of the regularized Green tensor consist 
of two parts: an oscillating and a decaying part, as a func- 
tion of |z — The decay occurs at a distance that is a 
tiny fraction of an optical wavelength. For A|z — zi| > 1, 
the regularized Green function approaches the unregu- 
larized one. If one would take the limit A — > cxd, then 
all the components in H24|l approach the unregularized 
components of Eq. (|18|l . and in particular the limit of 
the last term in Gg^ gives the delta function that made 
the regularization procedure necessary. However, A is 
kept finite for the moment, so that Gg^ has a finite term 
that grows with A. With this result, the Green-function 
regularization is complete and a theory of scattering by 
vector waves from plane scatterers can be set up. 



Here Gq^ stands for G^^ {k\\, Za, Za,uj), and similarly for 
the other components. The off-diagonal elements of the 
Green tensor are all zero when the position z is equal 
to Zq,, so that the equation can be solved for every com- 
ponent separately. By inserting this result into the LS 
equation for general z, one finds 



Eko-(k||,z,cj) = Eko-,o(k||, z,a;) 
-|-Go(k||,z,Za,w) • t(k||,w) • Eko-,o(k||,ZQ,a;), (26) 



where the T-matrix for scattering from a plane by arbi- 
trarily polarized light is given by 



t(k||,C.) 




















-y(a;)Gg- 



(27) 



The scattering of the s-polarization component of the 
light can be considered independently from the v and z 
directions, according to Eqs. (|26|) and (|27|) . It can be 
verified with Eqs. (HH EH) that since A > (w/c), the 
matrix component T*'' for all practical purposes is equal 
to the T-matrix for scalar waves, and the same holds for 
the Green tensor component Gff : the regularization was 
not necessary for s-polarized light and fortunately it does 
not affect the scattering properties of s-polarized light. 



C. T-matrix of a plane for vector waves 

The regularization entails that the Green function is re- 
placed by its regularized version in the LS equation (|21|l . 
For z ^ Za that equation becomes 



^k.(k|l 
^k.(k|| 



+V{uj) 




Eo 





























(25) 



The need for regularization did show up in the descrip- 
tion of scattering of p-polarized light, and there the cutoff 
might influence light scattering. Incoming p-polarized 
light is characterized by its amplitude i^o, wave vec- 
tor k, and its polarization state cr = f) = (fcz/fc)vk — 
(/c||/fc)z. Written out explicitly, the incoming field is 
EkCT,o(k|| , z, cj) = i?o(0, fcz/Zc, — A:||//c) exp(?fczz). For dis- 
tances far enough from the plane so that A|z — Zq| 3> 1, 
the term Gg^f in Eq. ^ falls off as A"! and Gg^f 
as exp(— A|||z — Zq|), so that for optical purposes these 
terms can be neglected. For finite very large A we arrive 
at the following effective description: 



E^{z)\ El[z)\ /Gr \ /T- 0\ E^iz^) \ 

E-{z) = E-oiz) \ + \ Gg'' GgM i?g(^a) , (28) 

E^{z) J \E^iz)J \0 Gg"Gg^y\ j\E^iz^)J 



7 



where the ss-component of the T-matrix is equal 
to V{uj) [1 — V{oj)Gq'']~^ , and analogously for the vv- 
component. The Green functions have arguments 
(k|| , z, Zq, cj). In this effective description, - where the 

T-matrix is denoted by T rather than T - the cutoff pa- 
rameter A does not occur anymore. The cutoff was ne- 
cessary in order to set up a scattering theory and it shows 
up in the elements of the scattering theory such as the 
T-matrix 1)2 7|l and the regularized Green function (|24|l . 
It does not show up in the electric field, and precisely this 
enables us to arrive at the effective description. Note also 
that the (large) value of Gq^ has become irrelevant. 

The effective description that is obtained here after a 
regularization is different from a theory where the delta 
function in Gq^ [see Eq. (|18d|l ] would simply be removed 
|35|. Leaving out the delta function in the LS equa- 
tion l|21|l will result in a nonzero T^^, in contrast with 
Eq. H28() . Furthermore, the T-matrix would have the un- 
wanted effect that the transmitted part of an incoming 
wave would not be parallel to the incoming wave. The 
conclusion is that a regularization of the Green function 
was necessary, even when in the end the cutoff could be 
sent to infinity. 

The equation H28(l defines a true mode of the electro- 
magnetic field in the presence of a single plane scatterer, 
in terms of a linearly-polarized incoming plane wave with 
arbitrary angle of incidence. This is not the complete set 
of modes. Other modes, not corresponding to an incom- 
will be discussed in Sec. 



mg wave, will be discussed m bee. (jlVBI 
single plane and for a crystal of planes. 



both for a 



D. Transmission and energy conservation 

The transmission of light through the plane can be 
found by choosing z > Za in Eq. H28|l . The transmitted 
wave can be expressed in terms of the incoming wave 
as Ekcr(k|| , z, w) — T(k||,a-') • Ekcr.o(k|| , z, w), with the 
transmission matrix 



Im T"^(k 



II' 



2 



'(k|h^)l^ 



(30) 



The most general T-matrix satisfying this requirement 
has the form 



T-(k||,c.) = -[F7i(k||,c.)-»/(2fc.)]" 



(31) 



where the optical potential _F's(k||,w) must be a real- 
valued function. 

For reflection and transmission of p-polarized light, 
only the matrix element T"" is important. Again, we 
are interested in the form that this matrix element can 
take when optical energy is conserved. This is the case 
when the z-component of the Poynting vector is the same 
before and after the plane. An incoming p-polarized 
plane wave gives the electric field (EHJ. With a MaxweU 
equation the accompanying magnetic field B can also 
be found. In Sl-units, and in terms of the complex 
fields E and B, the cycle-averaged Poynting vector is 
equal to Rc[E*(r,t) x B(r, t)]/(2^o) (3Zj- When a har- 
monic wave of frequency uj coming from z — ^oo scat- 
ters off the plane, the Poynting vector is proportional to 
1 - (A:^c/w)2|T™|V4 for (z < z„). At the other side of 
the plane one finds |1 — ikz{c/uj)^T^^ /2\^. By equating 
the two, the optical theorem for a plane that scatters 
p-polarized light is found to be 



ImT™(k||,c.) 



2{uj/cy 



|r™(k||,c.)|2 



(32) 



This differs from the optical theorem for s-polarized light. 
Also, the most general solution of the optical theorem is 
different: 



T™(k||,o.) 



ik. 



2{uj/cY 



(33) 



where the optical potential Fp(k||,a;) is real. 



r(k||,w) 



5(k||,w) 







r„t,(k||,u;) 
T2t,(k||,a;) 1 



(29) 



which has nonzero 



•30 



(k|h^) 



1 



elements 
1 

for i ~ s,v. Fur- 
thermore, Tzv{'k-\\,(^) — G^'"(fc||, Zq,, Zct,u;)T''"(fc||, w). 
Both for purely s-polarized and for purely p-polarized 
light, the transmitted electric field is a polarization- 
dependent scalar times the incoming electric field 
vector. 

Energy conservation puts a constraint (called 'optical 
theorem') on the form that the T-matrix of an elastic 
scatterer can take. The optical theorem for a plane that 
scatters scalar waves was found before ^5;] • Since s-waves 
map on scalar waves, the optical theorem for the ss- 
component of the T-matrix can be given immediately: 



E. T-matrix for N planes 

Now that the Green function and the T-matrix of a 
single plane are known, a multiple-scattering theory can 
be set up. Assume that there are N plane scatterers, 
placed at arbitrary positions. Assume them to be par- 
allel, so that s- and p-polarized light do not mix in the 
scattering process. 

In the general expression H12|) for the T-matrix of a 
complex dielectric in terms of its simple parts. Green 
functions are always sandwiched between T-matrices of 
scatterers at different positions. For unequal plane posi- 
tions Za and Zf3, the value of Go(k|| , z^, Zq,, a;) is finite 
and it can be taken to be equal to the unregularized 
Go(k|| , z^, Zq,, w), because different planes are at optical 
distances apart. Further regularizations are therefore not 
required in order to find the A^-plane T-matrix. 
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As shown in the Appendix, higher-order terms in the 
series (I12II correspond to higher-order matrix multiph- 
cations oi N x N matrices. The multiphcation prop- 
erty makes that for parallel planes, the series (|12|) can 
be summed exactly. In the Appendix it is shown how 
the summation can be done even in the general case that 
all planes may have different optical properties, and are 
placed at arbitrary non-coinciding positions. Here we 
specify that all planes are identical. This gives the central 
result of this section, the TV-plane T-matrix for scattering 
by vector waves: 

1 r ^ 

(34) 

Each a/3-component T^^"^ (k|| ,lu) is a 3 x 3 matrix; the 
only two nonzero spatial components are 

j..s,(JV) ^ T^qi^ _ >v^T-];^ (35a) 
j.vv.iN) ^ [Iw - W„T™];^ (35b) 

Here, I at is the N x N unit matrix. Arguments 
(k|| , Lu) were temporarily dropped for readability. The 
N'^ matrix elements {yVj)ap(k.\\,i^) are defined as (1 — 
'5q/3)Gq-' (k|| , Zq,, z^, w), for j — s,t). The calculation of 

T*-^-* (k|| , w) boils down to the inversion of an TV x 
matrix for the two transverse polarization directions se- 
parately. 

From now on, assume that the N planes are placed 
at regular distances from each other, with a spacing a 
between neighbors. The necessary matrix inversions in 
Eq. (|35l) can then be performed analytically for both po- 
larization directions. The s-wave case maps identically 
on the situation for scalar waves, for which the analytical 
inversion was discussed at length in 5]; for p- waves the 
inversion trick goes analogous and it will not be presented 
here. 

A result from the analytical inversion is that T-matrix 
elements and therefore the optical properties of the A^- 
plane crystal strongly depend on the Bloch wave vec- 
tors Ars(k||,ijj) and Arp(k||,a;). For p-polarized hght 
the Bloch wave vector is given by arccos(Cp)/a, with 
Cp = cos(fcza) -|- CpShi{kza); the constant Cp in terms 
of the single-plane T-matrix is 

ik^ [fc^|T''"|2 + 2(w/c)2lmr™] + 2fc3(w/c)2Rer™ 
fc2(Rer™)2 -I- [2(a;/c)2 fcJmT™]^ 

(36) 

In general, Cp is a complex constant. However, if the op- 
tical theorem H32(l holds, then the imaginary part of Cp 
becomes identically zero, and the single-plane T-matrix 
will be of the form (|33|l . Likewise, Kg is defined as 
arccos(Cs)/a for a quantity Cs that becomes real when 
the optical theorem Eq. (|30|l for s-polarized light holds 
In those cases, the expressions for Cs^p become rather 



simple, 

/F^(k||,tj)\ 

Cs = cos{kza)~l — — j sm{kza) (37a) 

, /'k^c^Fp{kn,uj)\ , , , 
Cp = cos(fc2a) - ( ^^-jf j ii\n{k^a). (37b) 

F. A model for the optical potential 

The most general T-matrices (I31|l and (I33II feature as 
yet unspecified optical potentials Fs_p. These should be 
real when energy is conserved, but for the rest they can be 
arbitrary functions with the frequency and the in-plane 
wave vector as variables. 

In 1^, plane scatterers were introduced as a simpli- 
fied model for dielectric slabs of finite thickness d and 
nondispersive dielectric function £{uj) = e. The optical 
potential for the plane scatterer in this model is obtained 
via the limiting process of making the thickness d of the 
dielectric slab smaller and increasing the polarizability 
(e — 1), while keeping their product constant and equal 
to the "effective thickness" Doff- (The quantity Dcs/a is 
called the "grating strength" in [I^jIIqI-) Following the 
same limiting procedure as in Isl , we find the optical po- 
tential Fs_p(k||,a;) = —V{lu) = _Deff(w/c)2, identical for 
the two polarizations. Spatial dispersion and anisotropy 
would have shown up in the optical potentials as a and 
k|| -dependence, respectively. These two phenomena were 
neglected already as early as in the wave equation 

In general, p-polarized light differs from s-polarizcd 
light in that the former will have a Brewster angle at 
which no light is reflected from a dielectric interface 
(ni — > n2). The Brewster angle 9b equals tan~^(n2/ni). 
In the limiting procedure for going from a finite slab- 
in-air to an infinitely thin plane-in-air, the dielectric 
contrast y/e/1 is going to infinity and consequently the 
Brewster angle becomes 90° in that limit. Therefore, in 
our limiting procedure, a plane scatterer will not have a 
Brewster angle at the same angle as the finite dielec- 
tric slab that one starts out with. In line with this, 
in 15] a single-plane reflection for p-polarized light was 
determined that is nonzero for all angles of incidence. 
The p-polarized propagating modes for a system of plane 
scatterers will therefore differ substantially from the cor- 
responding modes in a slab structure. The absence of a 
Brewster effect was also noticed in where the infinite- 
crystal version of the plane-scatterer model is treated. 



IV. OPTICAL MODES AND 
OMNIDIRECTIONAL MIRRORS 

A. Propagating modes 

The optical modes are the harmonic solutions of the 
wave equation (gj. With the solution of the T- 
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matrix, the modes that correspond to a nonzero incoming 
plane wave can be given explicitly as 

Ek^(k||,z,u;) = £;ocrke^'=^^ (38) 
+ Go(k|| ,z,z^,Lj)- tI'P (k|| , uj) ■ rr^Eae'^^'-^ . 

a,f3 

These propagating (or radiative) modes are labelled by 
the incoming wave vector k and polarization cri^. The 
sine of the angle of incidence (with respect to a vector 
normal to the planes) is equal to fc||c/w. The amplitudes 
of the s-polarized modes [with <Tk = (1,0, 0)] are identical 
to the corresponding amplitudes for scalar waves; the p- 
polarized modes [crk = (0, fc^/fc, — fc||/fc)] have no scalar 
analogues. 

For light with wave vector and frequency such that 
Cs,p > 1 [see Eq. (I37f) ]. the Bloch wave vector is purely 
imaginary for the elastic scatterers that we consider. 
Similarly, for Cs,p < —1, the Bloch wave vector equals 
TT plus an imaginary number. In both situations the light 
will feel a stop band, meaning that it will be 100% re- 
flected when falling on a semi-infinite system of planes. 
Otherwise, when — 1 < Cs.p < 1, the Bloch wave vector 
is real and light can propagate inside the crystal. More 
will be said about the Bloch wave vectors later in this 
section. 

Some plots of mode profiles will now be presented. As- 
sume that light comes in from the left. For perpendicu- 
larly incident light, there is no difference between s- and 
p-polarization. In Fig. ^ the mode profiles (or squared 
absolute values of mode functions) for s- and p-polarized 
light inside a ten-plane crystal are compared both for an 
incoming angle of 30° and for 60°. Fig.^a) shows that at 
an angle of 30° the mode profiles corresponding to both 
polarizations do not differ much yet. Both modes decay 
rapidly inside the crystal structure and are reflected (al- 
most) completely. The Bloch wave vectors are complex 
for both polarizations. Only for the s-wave the polar- 
ization directions of the incoming and the reflected wave 
are equal, so that the amplitude of its mode profile at 
the left side of the crystal is four times the amplitude of 
the incoming electric field. 

The situation is different at an incoming angle of 60°, 
as shown in Fig.^b): the mode profile of the s-polarized 
light again rapidly decays inside the crystal (and the cor- 
responding Bloch wave vector again has an imaginary 
part), whereas the p-polarized light can propagate inside 
the crystal and is transmitted almost completely (and 
the Bloch wave vector is real). For this frequency and 
incoming angle, the crystal is a good polarization filter. 

The mode profiles in Figs. ^a,b) of the s-polarized 
waves are continuous whereas p-polarizcd waves are dis- 
continuous at the positions of the planes. This reflects 
the boundary conditions: the tangential components of 
the electric fields must be continuous and the normal 
components must show a jump at a dielectric interface. 
The electric field of s-polarized light only has a tangen- 
tial component, while p-polarized light consists of both 



tangential and normal components. This explains the dif- 
ferences in the mode profiles for s- and p-waves. Notice 
that in our Green-function formalism, boundary condi- 
tions are automatically satisfied, whereas in related work 
based on transfer-matrix methods, boundary conditions 
must be considered explicitly [H, [la, [23 . 

Reflection as a function of frequency by the ten-plane 
Bragg mirror is plotted in Fig.[2|for both polarization di- 
rections. The reflection equals (1— |r|^), with |r|^ the 
(relative) transmitted light intensity. For light incident 
perpendicularly to the planes, both transverse polariza- 
tion vectors are equivalent and accordingly in Fig. [2Ia), 
the graphs for s- and p-polarized light overlap. Differ- 
ences between the two polarizations do appear for non- 
normal incidence. In Fig. the angle of incidence is 
60° . The red edges of the stop bands for s-polarized light 
move to slightly higher frequencies and the widths of the 
stop bands become larger. For p-polarized light the red 
edges of the stop bands shift to the blue much faster, and 
the faster so for larger angles of incidence. 

For scalar waves, a crystal of plane scatterers can be 
an omnidirectional mirror [^, which means that waves 
experience a stop bands for all angles of incidence. For 
vector waves, the crystal will only be an omnidirectional 
mirror if there are frequency intervals in which the crystal 
is an omnidirectional mirror both for s- and p-waves. 

As stated earlier, it is the Bloch wave vectors K that 
distinguish between light that can propagate inside the 
crystal (real K) and light that feels a stop band (when K 
is an imaginary number or tt plus an imaginary number). 
In our formalism, the Bloch wave vectors are the arc 
cosines of the constants Cg and Cp given in Eq. (|37|l . As 
is shown in detail in Q , these Bloch wave vectors show up 
in expressions for the A'"-plane T-matrix T^^' . It must be 
said that in the present T-matrix formalism it is not obvi- 
ous simply by looking at the equations that a stop band 
occurs whenever the Bloch wave vector has a nonzero 
imaginary part. Nevertheless, we conclude from our nu- 
merical calculations that the relation does exist. To give 
an example, of the four modes in the figures^a,b), only 
the p-polarized light incoming at 60° has a correspond- 
ing real Bloch wave vector. It is an interesting fact that 
the Bloch wave vector, the role of which is obvious in 
infinite crystals, also plays an important role in finite pe- 
riodic structures. This was already noticed before in the 
context of transfer matrix methods [53, [H, [1^ ; here we 
see the importance of the Bloch wave vector for finite 
periodic structures in a T-matrix formalism. 

For light of a frequency corresponding to a/ A = 0.5 
and planes with Lha = 0.46, the s-waves are reflected 
omnidirectionally [3 . In Fig. [3| we plot both constants 
Cs and Cp for this frequency, as a function of angle of inci- 
dence. Unlike for s-waves, for p-waves there are incident 
angles larger than the critical angle dc > 55° for which 
the values of Cp are between -1 and 1. Light incident 
with these large angles can propagate inside the crystal 
and therefore the crystal is not an omnidirectional mir- 
ror for this frequency. Actually, this information could 
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Position z (units of a) Position z (units of a) 



FIG. 1: Squares of absolute values of mode functions for s-polarized (solid lines) and p-polarized light (dashed lines), as a 
function of position. The light is scattered by a crystal of ten planes with D^h — 0.46a, separated by a distance a. Both modes 
correspond to light incoming from the left with a/A = 0.5. Figure (a): ^in — 30°; figure (b): 6'in = 60°. 




FIG. 2: Reflection off a ten-plane crystal, as a function of a/A, for s-polarized light (solid lines) and p-polarized light (dashed 
lines). Angles of incidence are 0° in (a) and 60° in (b). The planes have effective thickness Dcff ~ 0.46a and they are separated 
by a distance a. In (a), the graphs for s- and p-polarization overlap. 



already be read off from the mode profile of p-polarized 
light incident at 60° in Fig. ^b). The conclusion holds 
more generally: for larger D^s the critical angle 9c in- 
creases, but it can be shown by expanding Eq. Ij37b|) 
around ~ 90° that for every finite D^ff/a and a/A 
there always is a finite interval of angles corresponding 
to propagating p-polarized light. In conclusion, crystals 
of identical and equidistant plane scatterers can reflect 
vector waves in almost all (but not in all) directions. 



B. Guided modes 

Besides propagating modes there can also be bound 
modes that do not correspond to incoming light (see 
Sec. Bound modes can be found by solving the LS 
equation in the absence of an incoming field. In crys- 
tals of plane scatterers, bound modes are guided modes. 
They have imaginary wave vectors in the z-direction and 
they decay exponentially away from the planes. Their 
in-plane wave vectors fc|| are larger than lo/c. With each 



mode, be it of the propagating or guided type, a nonzero 
local density of states is associated. In the following, 
guided modes will be searched by looking for nonzero 
densities of states. (This is a method alternative to the 
one used in Q where guided modes of scalar waves were 
identified.) 

For vector waves, the local optical density of states 
N{r,Lu) is defined by 

-[(2c^)/(^c2)]ImTrG(r,r,c^), (39) 

so it is a scalar proportional to the trace over the imag- 
inary part of the Green tensor G(r,r,ti;). In planar ge- 
ometries the latter can best be found as an integral over 
the Green tensor in the plane representation: 

G(r,r,w) = J d2k||G(k||,^,z,w). (40) 

The local density of states can only be nonzero if the 
imaginary part of the integrand in (|40|) is nonzero. A 
guided mode manifests itself when diagonal elements 
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FIG. 3; Constants Cs and Cp as a function of angle of the 
incoming hght, for the parameters a/A = 0.5 and D^e — 
0.46a. Regions where — 1 < C < 1 correspond to propagating 
waves inside the crystal. 



of this integrand G(k|| , z, z, w) have nonzero imaginary 
parts for a certain fc|| > lu/c. For the crystals of plane 
scatterers the Green tensor directly follows from the 
Dyson-Schwinger equation: 



G(k||, z, z',cj) = Go(k||,z,z',cj) 



(41) 



N 



- X! ^Om, Za, UJ) ■Tfp\'ki\,Uj) ■ Com, Zf},z',L 
a,l3=l 



All three diagonal components of Go(k|| , z, z', ui) become 
real quantities for fc|| > w/c, and indeed there are no 
guided modes in free space. On the other hand, the 
off-diagonal elements Gg^ = Gg" become purely imag- 
inary when fc|| > Lu/c and the latter elements do show 
up in the diagonal elements of G. However, since they 
always show up in paired products, for example in the 
term Gg''T^^^'™G5^ they also give a real contribution 
to diagonal elements of G. The T-matrix elements are 
also real when fc|| > lu/c, except when the matrix has 
a pole. Therefore, all guided modes must correspond to 
poles of the s- or p-components of the A'^-planc T-matrix 



-(AT) 



Eq. (ESI] 



First the guided modes of a single plane will be deter- 
mined. There can be a guided mode when either T*'* 
or T"" in Eq. ^ has a pole. Now T'*" has a pole 
when 1 — F(ci;)GQ''(k|| , Zq, z^, w) vanishes. Using the same 
model for the optical potential as in Sec. 1111 11 we find 
the dispersion relation = Z)eff('^/c)^/2 for one and 
only one guided mode corresponding to s-polarized light. 
Here, k is the positive square root [fc^ — (w/c)^]^/^ for 
fell > Lo/c. A single-plane guided mode with this disper- 
sion relation was also found in 15], and in 5] for scalar 
waves. 

A pole of T™ occurs when 1 - V{uj)G^'' {k^, z^, z^^uj) 
vanishes, which is equivalent to the requirement that 
Dcsk/2 equals —1. Now in principle, D^s could be nega- 
tive when modelling a slab of negative dielectric function 



as a plane scatterer. However, in the physical situations 
that we are interested in, the effective thickness is real 

Therefore, there is no guided 



and positive (see Sec. lIII F|) 

mode corresponding to p-polarized light for a single plane 
scatterer. This result was also derived in where only 
the cases of a single plane and infinitely many planes were 
considered. 

At this point it is worthwhile to compare the guided 
modes of a dielectric slab (thickness d, dielectric constant 
e) in air with the guided modes of an infinitely thin plane 
with effective thickness DaS = (e— 1)^- For the slab, the 
number M of guided modes is the same for both polar- 
izations in the special case considered here, and equal to 

M 



1 



(42) 



Here, [X] stands for the largest integer smaller than or 
equal to AT. In the large- wavelength limit there is a single 
guided mode for each polarization direction. The second 
guided mode appears when a/A = a / {2\/ dDos) ■ For ex- 
ample, for d = 0.1a and Dcs — 0.46a, a second guided 
mode exists when a/A > 2.3. We are interested in fre- 
quencies around a/ A = 0.5, where the first stopband for 
normally incident light occurs (see Fig.[2Jl. For these fre- 
quencies, both the plane scatterer and the dielectric slab 
have a single s-polarized guided mode; the slab has a sin- 
gle p-polarized guided mode, whereas the plane scatterer 
has no such guided mode at all. 

Now we determine the guided modes of a finite crys- 
tal of N parallel and equidistant planes, using the same 
Green-function method as for the single plane. First look 
for the poles of the component T^^'''^-'(k|| , w). This is 
easy, because this component is identical to the A^-plane 
T-matrix r(^) for scalar waves, for which it was found 
that there are at most N guided modes in a crystal of N 
planes |^. For an infinite number of planes, the guided 
modes form a band, as was found in [l^ [T^ • Now look 
for guided modes corresponding to p-polarized light. The 
poles of the component T'^^'^^^ (k|| , lu) occur when the de- 
terminant det[(T''"'^^))~^] is equal to zero. An expres- 
sion for this determinant can be found just like was done 
for scalar waves in . The result is that for p-polarized 
light there are guided modes if the following equation has 
nontrivial solutions ui{k): 



sm[{N + l)Kpa] 



KF--2{uj/cf 
2(w/c)2 



e-"" sm{NKpa) = 0. 

(43) 
times 



The Bloch wave vector Kp is still defined as 
the arccosine of Cp, which reads Cp — cosh(Ka) -f 
[{nc^F)/{2uj'^)] sinh(Ka) in terms of k. 

Eq. H43|) should lead to the dispersion relations uj{k) 
for the guided modes, if they exist. When increasing 
the frequency, new guided modes would appear that at 
first are only just captured by the structure so that k = 
0^. It is therefore convenient to count the guided modes 
in the small-K limit. Let the constant x be defined as 
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a^/T^^^\?JJJlLu^. Then Cp can be written up to second 
order in k as 

Cp = l + (xk)V2 + o(k3) = cosh(xK) + o(k^). (44) 

To the same order in k, the Bloch wave vector Kp be- 
comes equal to ixf^/a- Therefore, solutions of Eq. 
will only exist when sinh[(iV + 1)xk] equals sinh(iVxK), 
or equivalently when % = 0. Since F is taken to be 
£'oft(<^/c)^, there are guided modes for p-polarized light 
whenever 1 + Dcs/a = 0. So in finite crystals of N 
planes with a positive effective thickness, there are no 
guided modes corresponding to p-polarized light. This is 
in agreement with the result iust obtained for the single 
plane, and with the results in [l5Ll2(]| for infinite crystals. 

In conclusion, there are at most N guided modes in 
the finite crystal of plane scatterers, all modes corre- 
sponding to s-polarized light. The comparison of the 
single plane and the single slab indicates that s-wavcs in 
plane scatterers are a good model for s-waves in slabs, at 
least for frequencies around the first stopband. For the 
p-polarized guided modes, the conclusion must be that in 
the finite slab structures there are guided modes which 
have no analogues in the crystal of plane scatterers. 



V. SPONTANEOUS EMISSION 
A. Application to layered dielectrics 

In free space, the spontaneous-emission rate Tq of an 
atom with dipole moment fi and transition frequency 
equals fx^^l^ /{Sirheoc^). When embedded in an inhomo- 
geneous dielectric, the rate F will in general be different, 
as can be found with Fermi's golden rule js^ 

F(/x, R, f]) = ^ V ^l/x • E;(R)|2 S{uji - n). (45) 

The E; are the normal-mode solutions with eigenfre- 
quencies uji of the wave equation The spontaneous- 
emission rate can alternatively be expressed in terms of 
the Green function of the medium 



Im[/x-G(R,R,17)-/x]. (46) 



(See ^3 for early derivations of this relation; in j4lj a 
modern derivation is given for inhomogeneous and ab- 
sorbing dielectric media.) In Eq. H4t)|l . G is the classi- 
cal dyadic Green function of the electric-field wave equa- 
tion Q). For homogeneous dielectrics, it is known that 
the total Green function is the sum of a transverse part 
that describes radiative decayand a longitudinal part de- 
scribing nonradiative decay HJ. Here, the Eqs. and 
(|46|l are equivalent, because nonradiative decay is absent 
for dielectrics with real dielectric functions. 

Layered dielectrics (not necessarily plane scatterers) 
are translation invariant in two directions, which can 



be chosen to be the {x, y) directions. Spontaneous- 
emission rates will only depend on the z-coordinate of 
the atomic position R — (x,y,z). It is then easiest to 
first calculate the Green function in the plane representa- 
tion G(k|| , z, z, il). This Green function must be Fourier 
transformed back to real space as in Eq. 14U|) in order to 
find the local Green function of Eq. (|46|l that determines 
spontaneous-emission rates. 

A slight complication in doing the integration H4U|) 
is that the plane representation for G(k|| , z, z, fi) is co- 
rotating with the incoming wave vector k|| , a variable 
that must now be integrated over. A fixed basis {x, y, z} 
is needed instead and it is chosen such that the atomic 
dipole becomes (/ia;,0,/iz) in the new representation. 
Write the two-dimensional integral / d^k|| in polar co- 
ordinates as jQ°°dfc||fc|| J^^ dk\\. After doing the angular 
integral, only diagonal elements of the dyadic Green func- 
tion survive. 

The total spontaneous-emission rate is the sum of two 
contributions, the perpendicular and the parallel decay 
rate 



(47a) 



T^{z,n) = -^^Im^°°dfc||fc|| ^^^ + ^V b) 

(Green-function arguments (fc||,z,z,i7) were again 
dropped.) The parallel decay rate has a contribution 
both from s- and p-polarized light (through ff** and G^"" , 
respectively) whereas the perpendicular decay rate only 
has a p-polarized decay channel (through G^^). Notice 
that the (real) delta- function term in Gq^ does not play 
a role in the emission rates. The spontaneous-emission 
rates in Eq. H47a|) and Il47b|l are integrals over all possible 
lengths of the in-plane wave vector. Both rates can be 
subdivided into a propagating-mode (or radiative-mode) 
rate corresponding to the integration of from to f2/c, 
and a guided-mode rate which is the integral from Jl/c 
to infinity. 



B. Spontaneous emission near plane scatterers 

The general expressions obtained in Sec. lV Al for spon- 
taneous emission in layered structures will now be ap- 
plied to crystals of plane scatterers. Combine the gen- 
eral expressions (|47|) for spontaneous-emission rates in 
layered dielectrics with the Green functions in the plane 
representation that were determined in Eq. H41|) for a 
crystal of plane scatterers. Because of the absence of 
p-polarized guided modes, the parallel decay rate near 
plane scatterers can be subdivided into three (instead of 
four) parts: an s-polarized radiative-mode rate (sr), a 
p-polarized radiative- mode rate (pr), and an s-polarized 
guided-mode rate (sg). Again, because of the absence of 
p-polarized guided modes, the perpendicular decay rate 
F^ is purely radiative. Here is a list of the nonzero partial 



13 



decay rates: 

rp^(z,f7) = 
r,(z,i7) = 



entations. For a dipole perpendicular to the planes, 



2n vTt)^'"/.?'''-''^"**' 



To be precise, in Eq. (|48d(l it was used that the ten- 
sor element G^^ has a vanishing imaginary part (leading 
to a vanishing contribution to the density of states) for 
/c|| > ri/c; for the same reason, there is no guided- mode 
rate analogous to Eq. H48c|) corresponding to G""" . These 
properties were found in Sec. IIV Bl 

With all partial emission rates spelled out now, we first 
study spontaneous-emission rates near a single plane, for 
which the Green function in the plane representation 
(|41|l features the single- plane T-matrix of Eq. H28I) . In 
Fig. ^a) , spontaneous-emission rates as a function of po- 
sition are plotted for DeS = 0.46a. For both orientations 
of the dipole, far away from the plane the rate approaches 
the free-space value. Close to the plane, F^; is larger than 
Fq, but it consists of a rate into propagating modes that 
is less than Fp, and a guided-mode rate. Close to the 
plane, F^ is larger than Fq, but the maximal values of 
(the purely radiative) F^ occur somewhat away from the 
plane. 

The contribution of radiative and guided s-waves for 
an atom with fj, = fx^ is the same as for scalar waves with 
"scalar dipole moment" /i, but since the total decay rate 
Fq for vector waves is larger than for scalar waves, the 
relative contributions of s-waves to F/Fq are smaller for 
vector waves (by a factor 3/4). 

In Fig.^fb), the same rates are plotted, this time for a 
plane with DcS = 10a that refiects light almost ideally: 
near the plane, F^ is almost twice Fq. The maximum 
values of F^ still occur away from the plane, although 
this has become invisible in Fig. 2Jb). The propagating- 
mode part of F^; has decreased and is practically zero on 
the plane. The partial emission rate into the guided mode 
has a much larger (but finite, not shown) amplitude near 
the plane. The other prominent difference in the two 
figures is that the 'spike' in the emission rates due to the 
guided modes has become much narrower. Indeed, from 
Eqs. (|39I41() . it follows that the guided-mode rate decays 
exponentially away from the plane like exp(— 2k^^^ jzl). 
It follows from the dispersion relation for k'^^'' that was 
obtained in Sec. IIV 51 that an increase in a/A or in Dcs 
will give narrower spikes. 

In the limit that the atomic position z becomes equal 
to the plane position Zq, = 0, the spontaneous-emission 
rates can be calculated analytically for both dipole ori- 



F,(z„,f7) =Fo 2 



3 \ 3 



1 + 



arctan(^) 



where the dimensionless parameter ^ is defined as 
7rZ)off/A. Similarly, for a dipole parallel to the plane, 
the three partial contributions to the decay rate can also 
be expressed in terms of the parameter ^ alone: 



rr(;^o,f^) 
K'{z^,n) 



-To [l-earctan(l/0] 



4^2 

37rC 



Fn 



1 - 



arctan(^) 



ro- 



(50a) 
(50b) 
(50c) 



In Fig. O the relative rates are plotted as a function of ^. 
The results can be checked in two limiting cases: if Dcs = 
there is no plane and then indeed both F^ and Fa; are 
equal to the free-space value Fq. The other limit is that 
of a perfect mirror, when Z^cff (and consequently ^) is 
sent to infinity. This limit is not visible in the figure, but 
the limiting values are Fz/Fq = 2 and Tx/Tq = 0. These 
values indeed agree with the well-known emission rates 
for atoms near perfect mirrors [2lll2^l23j |. Emission rates 
into guided modes vanish in the perfect-mirror limit; at 
z = Za this is only the case by sending DoS to infinity 
before putting z equal to Zq = 0. This completes the 
discussion of emission rates near a single plane. 

Now consider emission rates inside and near a crystal 
of a number N of plane scatterers. Results will be pre- 
sented for N ~ 10. In Fig.|H{a-d), orientation-dependent 
spontaneous-emission rates are plotted for several fre- 
quencies. For clarity in the pictures, the positions of the 
planes at a, 2a, 10a are not shown as vertical lines this 
time. The most striking difference between F^: and F^ 
is that Fa: becomes very spiky near the planes, because 
only parallel dipoles can couple to the s-polarized guided 
modes (and because p-polarized guided modes are absent 
in our model). 

As the frequency increases when going from Fig. Efa) 
to Fig. Eld), Fa: becomes more spiky, because the partial 
emission rate into the guided modes (the difference be- 
tween the solid and the dotted lines in Fig. |SJ) becomes 
more concentrated near the planes. The term 'concentra- 
tion' is appropriate here, because the maximum ampli- 
tudes near the planes are higher for narrower spikes. The 
same effect was observed for the single plane in Fig. 01 a- 
b), where the frequency is kept constant and Dcs is in- 
creased instead. 

The purely radiative-mode rate F^ on average increases 
due to the presence of the planes, whereas the radiative 
part of Fa: on average decreases. The same behavior oc- 
curs near a single plane in Fig.0|and for a perfect mirror. 

Figure n showed that the optical modes of p-polarized 
light have discontinuities at the plane positions. There 
are also discontinuities in the spontaneous-emission rates 
(not for a single plane, for symmetry reasons), but these 
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Position (units of a) Position (units of a) 



FIG. 4: Spontaneous-emission rates of dipoles near a single partially transmitting plane scatterer, relative to To. The wavelength 
is chosen such that a/\ — 0.5. Solid lines correspond to total spontaneous-emission rates Fx for dipoles parallel to the plane, 
dotted lines are radiative contributions to Tx, and the dashed lines denote Tz- In (a), the effective thickness DcS of the plane 
equals 0.46a and in (b) Deft ~ 10a. 




FIG. 5: Spontaneous-emission rates of dipoles at the position 
of the single plane scatterer, as a function of the dimensionless 
parameter ^ — TrDoff/A. The rate Tx (solid line) is the sum 
of a rate into radiative (dashed line) and into guided modes 
(dotted line). The purely radiative- mode rate F^ is the dash- 
dotted line. 



are too small to be visible in Fig. O It can be understood 
that they are small from the fact that the discontinuities 
per mode are averaged in the emission rate. 

The dotted lines in Fig. El are the radiative parts of 
Fa;. These are similar to the radiative- mode rates for 
scalar waves 0, but not identical, since in F^, not only 
s-polarized but also p-polarized light contributes, accord- 
ing to Eq. I)47b|l . In particular, far away from the planes, 
the emission rate of dipoles parallel to the planes con- 
sists of 75% s-polarized and 25% p-polarized light. (To 
be sure, light emitted by perpendicular dipoles is 100% 
p-polarized for all layered dielectrics, see Eq. (|47a|l .') 

For 'large enough' photonic crystals, one expects in- 
ner unit cells to have optical properties similar to unit 



cells in the infinite crystal. Are the ten-plane crystals 
large enough? This depends on the properties of a single 
plane. In two extreme cases, the crystal size does not 
matter: when the individual planes do not reflect any 
light in any direction (Dos/a = 0), or when they reflect 
all light {Dcff/a = oo), one finds the same emission rates 
in finite and in infinite crystals, because for Dcs/a — 
we have the free-space case and for Dgg/a — oo all unit 
cells are optically disconnected. Only in intermediate 
cases (0 < Des/a < oo) can finite photonic crystals 
have an appreciable and unit-cell dependent influence on 
spontaneous-emission rates. 

In the intermediate cases (we assumed Dcs/a = 0.46), 
the single-plane reflection also depends on the frequency 
of the light. The planes reflect light better at higher 
frequencies, because material dispersion was neglected. 
It appears that for the lowest frequency considered in 
Fig. El a/A = 0.2 in (a), the emission rates in inner unit 
cells already are influenced considerably by the crystal. 
The rates vary in neighboring inner unit cells, which indi- 
cates that the rates have not yet converged to the infinite- 
crystal values. For higher frequencies, say a/A = 0.6, 
individual planes reflect light much better. In the cor- 
responding Fig. jSJd), all inner unit cells look alike and 
emission rates have converged. 

Now consider a parallel dipole at a fixed position, very 
close to a plane in an inner unit cell. From Fig. El one 
can also appreciate how the guided modes will influence 
the (frequency-dependent) emission rate of a dipole very 
close to a plane. For frequencies a/A = 0.4 and higher. 
Figs. Efb-d) show that emission into propagating modes 
modes is negligible compared to emission into guided 
modes. The dipole falls inside the guided-mode spike 
near the plane. As the frequency increases, the maxi- 
mum amplitude of the spike increases. With the dipole 
still well inside the spike, the emission rate of the dipole 
will increase as well. When increasing the frequency fur- 
ther, the spike becomes so narrow that the dipole ends 
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FIG. 6: Spontaneous-emission rates (solid line) and Fz (dashed) near a ten-plane crystal. For all planes, D^s ~ 0.46a. The 
dotted line is the radiative part of Tx- The figures (a)-(d) correspond to four frequencies: (a) a/A — 0.2; (b) a/A — 0.4; (c) 
a/A = 0.5; (d) a/A = 0.6. 



up in one of the wings of the spike, until the dipole finds 
itself completely outside it. This will cause the emission 
rates into guided modes to drop at higher frequencies. 
The combined effect is a peak in the frequency-dependent 
emission rates. Indeed, in 19] the dipole emission rate (or 
more precisely, its s-wave component) for infinite crystals 
as a function of frequency shows a pronounced peak for 
dipole positions z near a plane (see Fig. 3 in |l3)- We 
can unambiguously attribute this peak to emission into 
the guided modes. 

We can also understand how an emission peak will de- 
pend on Doff and on the distance to the plane. We have 
seen that spikes are narrower, with higher amplitudes, for 
larger DcS and for higher frequencies. When assuming 
larger DcS, a dipole at fixed distance will feel a guided- 
mode enhanced emission rate for lower frequencies. How- 
ever, the dipole will also at lower frequencies begin to fall 
outside the range of the spike. This explains why for fixed 
dipole position and increased Dcff, the emission peak has 
a larger amplitude and attains its maximum at a lower 
frequency, precisely as seen in Fig. 3 in . Similar rea- 
soning suggests that for a parallel dipole a bit further 
away from a plane but still close to it (while keeping D^s 
fixed), an emission peak will occur at lower frequencies, 
with lower amplitude. 

What can be appreciated best in Fig. Efa) is that for 
parallel dipoles the total emission rate converges faster 



than the radiative and guided-mode partial rates sepa- 
rately. This can be related to a kind of mode competi- 
tion between s-polarized propagating and guided modes 
that we also found for scalar waves [Bj. On the other 
hand, for perpendicular dipoles all emission is into ra- 
diative modes, so mode competition is absent there and 
convergence sets in earlier. 

For scalar waves the ten-plane structure could act as 
an omnidirectional mirror, whereas in Sec. IIV Al it was 
found that it is not an omnidirectional mirror for vec- 
tor waves. Correspondingly, the radiative-mode LDOS 
for scalar waves at a/A = 0.5 dropped down to (almost) 
zero inside the ten-plane omnidirectional mirror, whereas 
the emission rates in Fig. ^c) show that the radiative 
LDOS for vector waves stays nonzero inside the crystal. 
In the inner unit cells, dipoles parallel to the planes emit 
predominantly guided light. The small amount of light 
that leaves the structure is strongly p-polarized. This 
is the case around a/A = 0.5 only, where s-polarized 
light is omnidirectionally reflected. Such strongly polar- 
ized emission is not a peculiarity of the plane-scatterer 
model, because it will also occur for a real Bragg mirror 
whenever light of only one of the two polarization direc- 
tions is omnidirectionally reflected. In the other plots 
in Fig. for higher and lower frequencies, the radiative- 
mode parts of are the sums of emission rates into both 
polarization directions. 
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VI. CONCLUSIONS, DISCUSSION, AND 
OUTLOOK 

A theory was set up for the multiple scattering of vec- 
tor waves by parallel planes, thereby generalizing pre- 
vious work for scalar waves to the more interesting but 
also more complicated case of light waves. Unlike for 
scalar waves, the Green function had to be regularized. 
This was accomplished by introducing a high-momentum 
cutoff. An effective scattering theory emerged with a 
nonzero T-matrix that no longer depends on the cutoff. 
The T-matrix and Green-function formalism turned out 
to be very convenient for the calculation of propagating 
and guided modes, as well as spontaneous-emission rates, 
of finite photonic crystals of plane scatterers. 

A non-absorbing plane scatterer satisfies a separate op- 
tical theorem for s- and p-polarized light. The radiative 
and guided modes of s-polarized light could be mapped 
onto modes for scalar waves. The s-polarized light has 
continuous modes, whereas p-polarized modes have dis- 
continuities at the plane positions. 

Throughout the paper, we have stressed the similarities 
and differences of optical properties of a plane scatterer 
or a crystal of plane scatterers as compared to the cor- 
responding dielectric slab structures, ft turns out that 
^»-polarized waves differ more in the two cases than s- 
polarized waves. Firstly, because propagating p-polarized 
modes are different in the two cases because the Brewster 
angle is 90° for plane scatterers. Secondly, p-polarized 
guided modes in finite slab structures have no analogues 
for plane scatterers. This was also found in for the 
single plane and in 0, |23| for the infinite crystal. 

Unlike for scalar waves, equidistant and identical plane 
scatterers can not be an omnidirectional mirror for all 
vector waves. Such omnidirectional mirrors consisting of 
dielectric layers do exist. For layered media, at least three 
refractive indices are required in order to prevent com- 
plete transmission of jo-polarized light in the Brewster- 
angle direction 3.. H. 11^ . 

Omnidirectional reflection is a property of a dielectric 
for external light sources. In this paper, omnidirectional 
reflection could be related to the emission properties of 
atomic light sources from within the finite crystals. The 
graphs of emission rates (Fig. ^ show that the emission 
by dipoles oriented parallel to the planes is affected much 
more strongly by the planes than emission by perpendicu- 
lar dipoles. This is a characteristic of the plane-scatterer 
model, because the absence of p-polarized guided modes 
is responsible for much of the difference. In the frequency 
interval where s-polarized light is omnidirectionally re- 
flected, all light that exits the crystal after a spontaneous- 
emission process will be p-polarized, irrespective of the 
orientation of the emitters. Still, the major fraction of 
the light will be emitted into guided modes and stay in- 
side the crystal. 

For low frequencies, the single-plane reflectivity is 
lower, emission rates are less affected by the crystal, and 
flnite-size effects are appreciable also in the inner unit 



cells of the ten-plane crystal. For higher frequencies, 
planes reflect light better and emission rates are more 
strongly modifled. In the inner unit cells, the emission 
rates converge faster to the values of the unit cell of an in- 
finite crystal. If the infinite crystal has larger variations 
in the emission rates inside a unit cell, then a smaller 
finite crystal is needed to converge to this result. 

We argued that the guided modes will give rise to a 
peak in the frequency-dependent emission rate of a par- 
allel dipole close to a plane. We also reasoned that the 
peak will shrink and shift to lower frequencies when ei- 
ther the effective thickness or the distance to the plane 
is increased. Indeed, the occurrence of a peak and its 
dependence on the effective thickness agree with numer- 
ical calculations on infinite crystals in p^ : some more 
numerical work is needed to corroborate our prediction 
of the distance dependence. 

What other purposes can plane scatterers serve in the 
future? The finite photonic crystals that can be made 
with them have one-dimensional periodicity only, yet 
light propagation in all three dimensions for all polar- 
ization directions is properly taken into account. The 
number of planes can be chosen at will and further advan- 
tages are that all optical modes and the complete Green 
tensor can be determined. 

Our calculations can be extended to situations where 
not al planes are identical; or one could allow light ab- 
sorption or gain in the planes by giving the effective 
thickness a complex value; the number of planes per unit 
cell could also be increased to more than one. Such cal- 
culations are possible in our formalism because the re- 
sults (|34|) and (|35|l for the T-matrix were generalized 
in Eq. lf5B|) of the Appendix to planes chosen at arbi- 
trary positions, each with a different T-matrix TQ,(k|| , lu), 
in other words with a different effective thickness. The 
model can therefore also be used to study the effects of 
disorder in the positions or in the optical properties of 
the planes on the spontaneous-emission rates of embed- 
ded atoms. 

It would also be interesting to study models of finite 
photonic crystals built up of non-parallel planes. This 
can be done, but numerically the model would become 
more involved, not so much because the planes have an 
overlap of measure zero, but rather because the plane 
representations for nonparallel planes will be different. 
Numerical calculations for nonparallel planes require the 
discretization of individual planes. 

Apart from extending it, one can extract other inter- 
esting observables from the model. For example, the 
knowledge of the complete Green function makes it pos- 
sible to calculate both far-field and near-field spectra of 
atoms embedded in the finite crystal. It would be in- 
teresting to study spectra near frequencies where the 
corresponding infinite crystal gives rise to a Van Hove 
singularity in the emission rates Crystals of plane 

scatterers can serve as a model environment to study the 
modification of several quantum optical processes when 
atoms are embedded in photonic crystals. Transient ef- 
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fects in the spontaneous-emission rates are but one exam- 
ple, thereby generahzing work done on a one-dimensional 
cavity formed by two planes 0| ■ Calculations are under- 
way that show photonic-crystal induced modifications of 
cooperative atomic processes. 



The third-order term of the expansion H12() is as Eq. H52|l . 



that is, the 



with Dai3 replaced by D^^ 
square of the NxN matrix D 
term features the cube of D, and so on. By summing all 
orders, we have 



7=1 —aj-'-'jIl, 

Similarly, the fourth-order 
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(54) 

where the TV-plane T-matrix T^^"^ (k|| , ui) is given by the 
series 



APPENDIX: DERIVATION OF N-PLANE 
T-MATRIX 
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D + D^ 
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(55) 



The goal of this appendix is to derive the expres- 
sions (|a4l35|l for the TV-plane T-matrix T^^\ by sum- 
ming the infinite series in Eq. H12(l , where the Ta are the 
T-matrices of single planes: 



Tq(w) 



(2^)2 



d'kii 



(k||,tj)(za,k|||. (51) 



The TQ(k|| , uj) are 3x3 tensors, defined in Eq. (j^HJ- First 
allow all planes to have different TQ,(k|| , w); furthermore, 
allow the parallel planes to have arbitrary (but all differ- 
ent) positions z^. 

The first-order term in the expansion (|12|l is simply 
the sum of the of Eq. I|51(l ; the second-order term has 
the form 



1 



(2^)2 



N 

/ dk|| |k||,z„)T„(k||;w)-[D(k||,w)]^^(k||,Z0| 



a,/3=l 



(52) 

where the dot denotes the inner product of 3 x 3 ten- 
sors. Furthermore, the property (k|| , z|Go(w)|k'|| , z') = 
(27r)^(5^(k|| — k'||)Go(k||,z,z',w) of the free-space Green 
function was used. The matrix elements D^^ of the ma- 
trix D are 3x3 tensors which are defined as 



[D(k|h^)]„« = (l-'5a/3)Go(k||,z, 



Here, 
the 3 



I N is the NxN unit matrix and, 
X 3 unit tensor, so that [I (g) Ijv] 



as before, I is 

a/3 = Now 

in the plane representation, the Green tensor has the 
form as derived in Sec. IIII Al and the single-plane T- 
matrices Tct(k||,w) have the simple form as defined in 
Eq. with only T^' and T^^ nonzero. It follows that 
Dap, Da0, and all higher powers of D have only three 
nonzero spatial components, namely the ss-, vv-, 
the zw-components. As a consequence, T^^''(k|| , 
only two nonzero spatial components for all plane indices 
a and /3, namely the ss- and ww-component: 



and 
has 



T^^Vk,, 

' a/3 1*^11 



T^^il^-D^-; 

^) = [1^ - D™];i 

V 

(56) 

This generalizes the result of Eq. (|28|l for the single-plane 
T-matrices T^, where the same two components were 
found to be nonzero. The generalization is from one plane 
to N parallel planes placed at arbitrary positions, each 
plane possibly with a different effective thickness. In the 
special case considered in the main text that all planes 

T/3(k||,w). 

are identical, in other words if all Tct(k||,Ci;) are equal, 
(53) then Eq. simplifies into Eq. 
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